Probing NMSSM Scenarios with Minimal Fine- Tuning by Searching for Decays of the 

Upsilon to a Light CP-Odd Higgs Boson 
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Completely natural electroweak symmetry breaking is easily achieved in supersymmetric models if 
there is a SM-like Higgs boson, h, with mi, <J 100 GeV. In the minimal supersymmetric model, such 
an h decays mainly to bb and is ruled out by LEP constraints. However, if the MSSM Higgs sector is 
expanded so that h decays mainly to still lighter Higgs bosons, e.g. h — > aa, with Br(/i — > aa) > 0.7, 
and if m a < 2mb, then the LEP constraints are satisfied. In this letter, we show that in the next-to- 
minimal supersymmetric model the above h and a properties (for the lightest CP-even and CP-odd 
Higgs bosons, respectively) imply a lower bound on Br(T — > 7a) that dedicated runs at present 
(and future) B factories can explore. 



Low energy supersymmetry remains one of the most 
attractive solutions to the naturalness / hierarchy prob- 
lem of the Standard Model (SM). However, the mini- 
mal supersymmetric model (MSSM), containing exactly 
two Higgs doublets, suffers from the "/i problem" and re- 
quires rather special parameter choices in order that the 
light Higgs mass is above LEP limits without electroweak 
symmetry breaking being "fine-tuned" , i. e. highly sen- 
sitive to supersymmetry-breaking parameters chosen at 
the grand-unification scale. Both problems are easily 
solved by adding Higgs (super) fields to the MSSM. For 
generic SUSY parameters well-below the TeV scale, fine- 
tuning is absent [H and a SM-like h is predicted with 
irih Si 100 GeV. Such an h can avoid LEP limits on the 
tightly constrained e + e~ — > Z+b's channel if Br(/i — > bb) 
is small by virtue of large Br(ft, — > aa), where a is a new 
light (typically CP-odd) Higgs boson, and m a < 2m(, so 
that a — > bb is forbidden Q . The perfect place to search 
for such an a is in Upsilon decays, T — > 7a. The simplest 
MSSM extension, the next-to-minimal supersymmetric 
model (NMSSM), naturally predicts that the lightest h 
and a, hi and 01, have all the right features [E 0, @| • 
In this letter, we show that large Br(/ii — > aiai) implies, 



at fixed 



a lower bound on Br(T — > 701) (from now 



on, T is the 15* resonance unless otherwise stated) that is 
typically within reach of present and future B factories. 

In the NMSSM, a light ai with substantial Br(/ii -> 
a±ai) is a very natural possibility for m^-scale soft pa- 
rameters developed by renormalization group running 
starting from U(1)r symmetric GUT-scale soft param- 
eters Q. (See also d, 0] for discussions of the light ai 
scenario.) The fine-tuning-prcferred mh 1 ~ 100 GeV 
(for tan/3 > few) gives perfect consistency with preci- 
sion electroweak data and the reduced Br(/ii — » bb) ~ 
0.09 — 0.15 explains the ~ 2.3a excess at LEP in the 
Zbb channel at M b r ~ 100 GeV. The motivation for this 
scenario is thus very strong. 

Hadron collider probes of the NMSSM Higgs sector 
are problematical. The hi — > aiai — > 4r (2m T < 



m ai < 2m})) or 4 jets {m ai < 2m T ) signal is a very 
difficult one at the Tevatron and very possibly at the 
LHC [1, 0, [Hi El ■ Higgs discovery or, at the very least, 
certification of a marginal LHC Higgs signal will require 
a linear e + e~ collider (ILC). Direct production and de- 
tection of the ai may be impossible at both the LHC and 
ILC because it is rather singlet in nature. In this letter, 
we show that by increasing sensitivity to Br(T — > 701) by 
one to three orders of magnitude (the exact requirement 
depends on m ai and tan (3), there is a good chance of de- 
tecting the a±. This constitutes a significant opportunity 
for current B factories and a major motivation for new 
super- B factories. Even with ILC hi — > aiai data, mea- 
surement of Br(T — > 701) and ai decays would provide 
extremely valuable complementary information. 

As compared to the three independent parameters 
needed in the MSSM context (often chosen as /i, tan/3 
and Ma), the Higgs sector of the NMSSM is described 
by the six parameters 

A , K , A x , A K , tan/3 , /x eff , (1) 
where [i c ff = \{S) = As is the effective /x-term gen- 
erated from the XSH u Hd part of the superpotential, 
\A\SH u Hd is the associated soft-SUSY-breaking scalar 
potential component, and n and nA K appear in the ^ kS 3 
and ^ kA k S 3 terms in the superpotential and associated 
soft-supcrsymmetry-breaking potential. In addition, val- 
ues must be input for the soft SUSY-brcaking masses that 
contribute to the radiative corrections in the Higgs sec- 
tor and to the Higgs decay widths. Our computations for 
branching ratios and so forth employ NMHDEC AY [l2| ■ 
An important ingredient for the results of this paper is 
the non-singlet fraction of the ai defined by cos 9a in 



ai = cos 6 aA 



MSSM 



+ sin 9 a A s , 



(2) 



where As is the CP-odd Higgs boson contained in the 
unmixed S complex scalar field. The coupling of a\ to 
t + t~ and bb is then oc tan (3 cos 9a ', cos 9a itself has some 
tan/3 dependence with the net result that tan (3 cos 9a 
increases modestly with increasing tan/3. 
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In [3, i, i], we scanned over the NMSSM parameter 
space holding tan/3 and the gaugino masses M\^,^{jnz) 
fixed, searching for choices that minimized a numerical 
measure, F, of EWSB fine-tuning, i.e. of how precisely 
the GUT-scale soft-SUSY-breaking parameters must be 
chosen to obtain the observed value of mz after RG evo- 
lution. For F < 15, fine-tuning is no worse than 7%, 
and we regard this as equivalent to absence of significant 
fine-tuning. For the sample values of tan/3 = 10 and 
M h2 ,3 = 100, 200, 300 GeV (F only depends significantly 
on M3), to achieve the lowest F values (F ~ 5 — 6), the hi 
must be fairly SM-like and ~ 100 GeV is required; 
this is only consistent with LEP constraints for scenar- 
ios in which Br(/ii — > a±ai) is large and m 0l < 2mb- 1 
Crucially, for these scenarios one finds a lower bound on 
cos8a\, e.g. | cos 0a | ^ 0.04 at tan/3 = 10. As described 
in [B[, this is required in order that Br(/ii — » a\ai) > 0.7 
when m ai < 2mf,. 2 

Aside from EWSB fine-tuning, there is a question of 
whether fine-tuning is needed to achieve large Br(/ii — > 
a\a±) and m ai < 2mt> when F < 15. This was discussed 
in [B[ . The level of such fine-tuning is determined mostly 
by whether A\ and A K need to be fine-tuned. (For given s 
and tan/3, Br(/ii — > aiai) and m ai depend significantly 
only on A, k, ^4 a and A K ; all other SUSY parameters 
have only a tiny influence.) Since specific soft-SUSY- 
breaking scenarios can evade the issue of tuning A K and 
A\ altogether, in this study we do not impose a limit 
on the measures of A\, A K fine-tuning discussed in [B|. 
However, it is worth noting that we find that A\, A K fine- 
tuning can easily be avoided if m 0l > 2m T and cos 9 a is 
small and negative, e.g. near cosf?^ ~ —0.05 if tan/3 = 
10. In some models, the simplest measures of A\, A K fine- 
tuning are much larger away from the preferred cos 9a 
region and / or at substantially lower m ai values. 

We now turn to T — » 701. We have computed the 
branching ratio for this decay based on Eqs. (3.54), (3.58) 
and (3.60) of [l]| (which gives all appropriate references). 
Eq. (3.54) gives the result based on the non-relativistic 
quarkonium model; Eqs. (3.58) and (3.60) give the proce- 
dures for including QCD corrections and relativistic cor- 
rections, respectively. Both cause significant suppression 
with respect to the non-relativistic quarkonium result. In 
addition, there arc bound state corrections. These give a 
modest enhancement, rising from a small percentage at 
small m ai to about 20% at m ai = 9.2 GeV (see the refer- 



1 We should note that the precise location of the minimum in F 
shifts slightly as tan /3 is varied. For example, at tan /3 = 3 
(tan/3 = 50) the minimum is at roughly 92 GeV (102 GeV). 
However, for these cases the minimum value of F is only very 
modestly higher at m,f ll ~ 100 GeV, the LEP excess location. 

2 Also, as one approaches the U(1)r, A k , A\ — > symmetry limit, 
large Rr(h\ — » a\a\) is not possible. 
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o tan/S=10, ^=150 GeV, M 1A3 = 100,200,300 GeV 
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FIG. 1: Br(T -> 701) for NMSSM scenarios with various 
ranges for m ai : dark grey (blue) = m ai < 2m T ; medium 
grey (red) = 2m T < m ai < 7.5 GeV; light grey (green) = 
7.5 GeV < m ai < 8.8 GeV; and black = 8.8 GeV < m ai < 
9.2 GeV. The plots are for tan/3 = 10 and Mi A s(m z ) = 
100, 200, 300 GeV. The left plot comes from the A X ,A K scan 
described in the text, holding fj, e ft(mz) = 150 GeV fixed. 
The right plot shows results for F < 15 scenarios with m ai < 
9.2 GeV found in a general scan over all NMSSM parameters 
holding tan/3 and Mi, 2, 3 fixed as stated. 



ences in 



13(). 3 Form ai e [m Vb -2r^ 



»)i.+ 2r »?bL where 
50 MeV and T Vb - 50 MeV, the a\ mixes 
significantly with the 77b, giving rise to a huge enhance- 
ment of Br(T — > 701). We have chosen not to plot results 
for m ai > 9.2 GeV since we think that the old theoretical 
results in this region require further refinement. In Fig.[TJ 
we present results for Br(T — » 701) that are consistent 
with existing experimental limits 4 in two cases: (a) using 
a scan over A\,A K values holding /z e ff(mz) = 150 GeV 
and Mi ]2i3 (m 2 ) = 100,200,300 GeV fixed (in this scan, 
identical to that described in Ref. A and k are also 
scanned over and all other SUSY-breaking parameters 
are fixed at 300 GeV - results are insensitive to this 
choice and, therefore, representative of the whole pa- 
rameter space); (b) for the F < 15 points found in the 
NMSSM parameter scan described earlier. In both cases, 



In contrast, for a scalar Higgs, bound state corrections give a 
very large suppression at higher Higg s masses near Mr- 
We impose the limits of Fig. 3 of U3 . Fig. 4 of [l5|, and Fig. 7b 
of [id . The first two limit Br(T — > 7X), where X is any visible 
state. The first provides the only strong constraint on the m ai < 
2m T region. The third gives limits on Br(T — > 7X)Br(V — > 
t+t~) that eliminate 2m T < m ai < 8.8 GeV points with too 
high Br(T -> 701) (for m ai > 2m r , Br(oi -> t+t~) ~ 0.9). 
Since the inclusive photon spectrum from T decays falls as Ej 
increases, the strongest constraints are obtained for small m ai . 
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^=150 GeV, M 1AS = 100,200,300 GeV 
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FIG. 2: We plot Br(T — » 701) as a function of cos 9a for 
the A X ,A K scan, taking Mi,2, 3 (m z ) = 100,200,300 GeV, 
fi e a(mz) = 150 GeV with tan/? = 3 (left) and tan/3 = 50 
(right). The point notation is as in Fig. \T\ 



all points plotted pass all NMHDECAY constraints — all 
points have ~ 100 GeV, but avoid LEP constraints 
by virtue of Br(/ii — > a\a\) > 0.7 and m ai < Im^. 
For both plots, we divide results into four m ai regions: 



m 0l < 2m T , 2?n T < m ttI < 7.5 GeV, 7.5 GeV < 
m ttl < 8.8 GeV and 8.8 GeV < m ttI < 9.2 GeV. Fig. Q] 
makes clear that Br(T — > 701) is mainly controlled by 
the non-singlet fraction of the ai and by m ai . The 
only difference between the (a) and (b) plots is that 
F < 15 restricts the range of cos 9a to smaller magni- 
tudes (implying smaller Br(T — > 701)) and narrows the 
m ai bands. As seen in the figure, the cos 6*^ ~ —0.05, 
m Q1 > 2m T scenarios that can have no A\, A K tuning 
have Br(T — > 701) < few x 10 -5 . For general cos Oa and 
m ai , values of Br(T — ► 701) up to ~ 10~ 3 (5 x 10~ 3 ) 
are possible for F < 15 points (in the general A\,A K 
scan). In Fig. [3 points with Br(T — > 701) > few x 10~ 4 
(depending on m ai ) are not present, having been elimi- 
nated by 90% CL limits from existing experiments. The 
surviving points with m ai < 9.2 GeV can be mostly 
probed if future running, upgrades and facilities are de- 
signed so that Br(T — > jai) ~ 10 -7 can be probed. As 
stated earlier, predictions at higher m ai arc rather un- 
certain, but obviously Br(T — > 701) — > for m ai —>- My. 
To access higher m ai (but m ai < 2mb), T(2S) — > 701 
and Y(3S) — > 701 can be employed; computation of the 
branching ratios requires careful attention to 01 — yyj, mix- 
ing, which can lead to even larger branching ratios than 
for the T if m ai ~ m Vb . 

Results from the A\,A K scan with /z e g = 150 GeV 
and Mi j2 ,3 = 100, 200, 300 GeV are given in the cases 
of tan/3 = 3 and tan/3 = 50 in Fig. Note that al- 



most all tan/3 = 3 points that pass NMHDECAY and 
LEP constraints are consistent with existing limits on 
Br(T -»■ yax). To probe the full set of m ai < 9.2 GeV 
points shown, sensitivity to Br(T — > 701) < fewxl0~ 8 is 
needed. Conversely, for tan /3 = 50 a lot of the scan points 
consistent with NMHDECAY and LEP constraints are 
already absent because of existing limits and one need 
only probe down to Br(T — > 701) ~ 10~ 6 to cover the 
m 0l < 9.2 GeV points. 

We note that the points with small negative cos 9a 
(e.g. cos6U ~ —0.5 for tan/3 = 10) that are most likely 
to escape A\ j A K tuning issues are well below the existing 
limits from [3, EH, 0] for all m ai values for all three 



tan /3 choices. 
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However, none of the above analyses 
have been repeated with the larger data sets 
available from CLEO-III, BaBar, or Belle. Presumably, 
much stronger constraints than those we included can be 
obtained. Or perhaps a 7«i signal will be found. 

We expect that the best way to search for the NMSSM 
light ai is to use its exclusive decay modes, as this reduces 
backgrounds, especially those important when the pho- 
ton is soft. For m ai > 3.6 GeV and tan/3 > 1, the dom- 
inant decay mode is a\ — > t+t~. For example, Ref. [l9[ 
has proposed looking for non-universality in T — > 7t + t~ 
vs. T — > 7e + e~,7^ + /i~ decays. This would fit nicely 
with the \ow-F scenarios. For m ai < 2m T ,2m c the de- 
cay mode ai — * gg is generally in the range 20% — 30%, 
giving a contribution to T — > 753 at the 10 -4 -10 -6 level; 
the ss mode is typically larger. 

In the 7t + t~ final state, the direct jt + t~ produc- 
tion cross section is 61 pb. Using signal=background as 
the criterion, this becomes the limiting factor for branch- 
ing ratios below the 4 x 10 -5 level when running on the 
T(15), and below the 2 x 10~ 4 level when running on the 
T(35). To improve upon the latter, one can select a sam- 
ple of known T(IS') events by looking for dipion transi- 
tions from the higher resonances. The dipion transition 
gives a strong kinematic constraint on the mass differ- 
ence between the two T's. When running on the T(35), 
the effective cross section in T(35) — > 7r + 7r~T(lS') 
is 179 pb Q. 6 To limit Br(T -> 7 ai) < 10~ 6 , 
5.6 fb _1 /e would need to be collected on the T(3S), 
where e is the experimental efficiency for isolating the 



5 For a CP-odd a that decays into non-interacting states, there are 
further constraints available from Crystal Ball and CLEO Il7|l ; 
these only apply to the scenarios considered here if Mi is reduced 
to a very small value (as possible without affecting EWSB fine 
tuning) so that ai — > x?X? decays are significant. For example, 
at tan ft = 10, our low fine-tuning scenarios with Mi decreased to 
3 GeV can yield m~ g 2 GeV and Br(ai -> x°X°) 6 [0.15, 0.35]. 

x i 

(Generic scenarios with substantial Br(T — > 7ai)Br(ai — > x?Xi) 
were considered in |lSj.) 

6 This can also be done on the T(25) but the pions are softer, 
implying much lower efficiency. On the T(45) this transition 
has a very small branching ratio ^ 10 -4 . 
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relevant events. This analysis can also be done on the 
T(45), where the T(3S) is produced via ISR. The effec- 
tive ji SR T(3S) — ► 7/5ij7r + 7r~T(lS') cross section is 0.78 
fb. To limit Br(T 701) < 10~ 6 , 1.3 ab _1 / e would need 
to be collected. These integrated luminosities needed to 
probe Br(T — » 70.1) ~ 10 -6 would appear to be within 
reach at existing facilities and would allow discovery of 
the a\ for many of the favored NMSSM scenarios. 

Are there other modes that would allow direct ai de- 
tection? 
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Reference 

with a\ — > 77. This works if the a± is very singlet, in 
which case Br(ai — > 77) can be large. However, see 
[H and earlier discussion, a minimum value of | cos 9 a | 
(e.g. I cos I > 0.04 if tan/3 = 10) is required in or- 
der that Br(/ii — > (Xiai) > 0.7 and m 01 < 2m(,. For 
the general A\,A K scans with Bi(hi — > aidi) > 0.7 and 
m ai < 2m& imposed, Br(ai — > 77) < 4 x 10~ 4 with 
values near /eu> x 10~ 5 being very common. It is con- 
ceivable that a super- B factory could detect a signal for 
T — > 7a 1 — > 777 which would provide a very interesting 
check on the consistency of the model. 

Flavor changing decays based on b — ► sai or s — > da\, 
in particular B — > A s ai, have been examined in Q. All 
penguin diagrams containing SM particles give contribu- 
tions to the b — > sai amplitude that are suppressed by 
cos 9a / tan /3 or cos 9a/ tan 3 /3 (since up- type quarks cou- 
ple to the Amssm with a factor of 1/ tan/3). Ref. 
identifies two diagrams involving loops containing up- 
type squarks and charginos that give b — ► sai ampli- 
tudes that are proportional to cos 9a tan /3. However, 
the sum of these diagrams vanishes in the super GIM 
limit (e.g. equal up- type squark masses), yielding a tiny 
B — > X s a\ transition rate. Away from this limit, results 
are highly model-dependent. In contrast, the predictions 
for T — > 7«i depend essentially only on cos 9a, tan/3 and 
m ai , ah of which are fairly constrained for the low-fine- 
tuning NMSSM scenarios. 

If m ai < 2m c , J/tp — > 7<2i decay will be possible. How- 
ever, Br(J/ip — > 701) is ~ 10~ 9 (~ 10~ 7 ) for the smallest 
(largest) | cos 9a\ values in the standard A\, A K scan for 
tan/3 = 10, increasing modestly as tan/3 increases. 

Before concluding, we note that a light, not-too-singlct 
a\ could allow consistency with the observed amount of 
dark-matter if the Xi is largely bino and 2m~ ~ m ai . 

This is explored in [l^|. We found that these scenarios 
could provide a consistent description of the dark matter 
relic density in the case of a very light Xi- We report here 
that this can be coincident with the F < 15 scenarios (as 
well as the small negative cos 9a, m ai > 2m T scenarios 
that are the most likely to have small A\, A K fine-tuning). 
All that is required relative to the Mi = 100 GeV choice 
made for our scans is to decrease M\ to bring down m~ 

near \ m ai . Mi is an independent parameter that has 
essentially no influence on the value of the fine-tuning 
measure F so long as Mi < M3. 



In summary, aside from discovering the a\ in h\ — > 
aiai decays, something that will almost certainly have 
to await LHC data and, because of the unusual final 
state, might not even be seen until the ILC, it seems that 
the most promising near-term possibility for testing the 
NMSSM scenarios for which EWSB fine-tuning is absent, 
or more generally any scenario with large Br(/ii — > aiai) 



and 



< 2mb, is to employ the T — > 701 decay at 



either existing B factories or future factories. 
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